function C = C_plane_strain(E,nu)
C = E/((1+nu)*(1-2*nu))*...
    [   1-nu nu 0;
        nu 1-nu 0;
        0 0 (1-2*nu)/2    ];
end
